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Abstract 

We  illustrate  a  computational  approach  to  practical 
nonlinear  balancing  via  the  forced  damped  pendulum 
example. 


1  Introduction 

The  theory  of  balancing,  introduced  by  Moore  [1],  was 
extended  by  Scherpen  [2]  to  a  class  of  stable  finite¬ 
dimensional  systems  in  the  nonlinear  setting 

x  =  f(x)+g(x)u  (1) 

y  =  Hx)  (2) 

where  x  =  (xi,...,xn)  are  local  coordinates  for  a 
smooth  state  manifold,  /,  g,  and  h  are  of  class  C°°, 
/( 0)  =  0,  and  h( 0)  =  0. 

The  existence  of  efficient  computational  methods  has 
been  crucial  in  making  balancing  an  effective  and  im¬ 
portant  tool  for  model  reduction  of  linear  systems.  The 
balancing  change  of  coordinates  can  be  computed  us¬ 
ing  established  matrix  equation  and  decomposition  al¬ 
gorithms.  In  contrast,  Scherpen’s  procedure  for  non¬ 
linear  balancing  presents  computational  difficulties.  In 
this  paper  we  address  these  difficulties  and  illustrate 
computational  methods  by  means  of  an  example. 


2  Balancing  for  the  Forced  Damped  Pendulum 

We  model  a  damped  pendulum  with  torque  input  and 
measured  angular  position  by 

X\  =  X2  ;  X2  =  —  asin(xi)  —  bx2  +  c u  (3) 
V  =  x  i  (4) 

where  we  choose  a  =  3,  b  =  1,  and  c  =  0.1. 
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The  main  objects  involved  in  Scherpen’s  theory  are  the 
controllability  and  observability  energy  functions,  Lc  : 
R"  — >  R  and  La  :  Rn  — »  R,  respectively,  defined  by 

1  r° 

Lc(x0)=  min  -  /  \[u{t)\\2  dt  (5) 

U  £E  2  ( —  OO,  0)  2  J- oo 
x(—oo)  =  0 
x(0)  =  x0 

and 

1  f°° 

Lo(x 0)  =  -  j  Ut)\\2dt,  (6) 

x(O)  =  xq  ,  u(t)  =  0 ,  0  <  t  <  oo. 

Computation  of  Lc  via  (5)  requires  solution  of  an  opti¬ 
mal  control  problem  at  each  point  on  a  state  space  grid. 
One  alternative  is  to  numerically  solve  the  Hamilton- 
Jacobi  type  partial  differential  equation  (PDE)  given 
by  [2]  eqn.  (12).  Both  of  these  methods  are  computa¬ 
tionally  intensive  (if  a  numerical  solution  can  be  gener¬ 
ated  at  all)  and  require  apriori  knowledge  of  /  and  g. 
We  offer  an  empirical  approach  which  could  be  used  in 
the  absence  of  accurate  models. 

In  the  linear  case,  Lc{x)  =  \  xT  Wc-1  x,  is  quadratic 
with  Wc  denoting  the  controllability  Gramian  matrix. 
Injection  of  Gaussian  white  noise  at  the  input  termi¬ 
nals  produces  a  Gaussian  asymptotic  probability  den¬ 
sity  function  (PDF)  for  the  state 

Poo(x )  =  [(27r)"det(f?oc,)]^1/2exp(-Tc(a:))  (7) 

where  I?oo  =  Wc  denotes  the  asymptotic  state  covari¬ 
ance.  Therefore,  in  the  nonlinear  setting,  we  estimate 
Lc  via 

Lc(x)  =  — log(p00(a;))  +  constant.  (8) 

This  suggests  an  empirical  Monte-Carlo  approach  in 
which  the  system  response  to  Gaussian  white  noise  is 
recorded  for  a  large  number  of  sample  paths,  and  the 
data  is  histogrammed  on  an  appropriate  grid  to  com¬ 
pute  a  PDF  and  estimate  for  Lc.  The  empirical  ap¬ 
proach  is  also  used  for  determination  of  La  via  numeri¬ 
cal  integration  of  (6)  where  the  input  signal  is  null  and 


the  initial  condition  Xo  takes  each  value  in  the  state 
space  grid.  Figures  1-3  show  results  for  the  pendulum 
system  (3-4). 


Figure  1:  Probability  density  function  for  asymptotic 
state  of  pendulum  system. 


Figure  2:  Estimated  controllability  function  for  pendulum 
system. 


Figure  3:  Observability  function  for  pendulum  system. 

The  existence  of  the  balancing  transformation  is  proved 
using  the  Morse-Palais  lemma  (see  [3,  4])  which  as¬ 
serts  that  there  exists  a  change  of  coordinates  under 
which  a  function  is  locally  quadratic  in  a  neighbor¬ 
hood  U  of  a  nondegenerate  critical  point  at  0,  i.e., 
/( x)  =<  A(j){x ),  (f>(x)  >,  with  (f)  a  diffeomorphism 
defined  on  U.  Computation  of  (f>  requires  the  decom¬ 
position  f(x)  =<  H( x)x,  x  >  (see  [3]).  It  is  shown 
in  [4]  that  <j>(x)  =  C{x)  x  where  C(x)2  =  B( x)  = 
H( 0)-1  H{x).  Such  a  C  exists  on  a  neighborhood  of 
0  since  B(0)  is  the  identity  and  a  square  root  function 
is  defined  in  a  neighborhood  of  the  identity  operator  II 
by  a  convergent  power  series.  The  algorithm  we  use  is 

So  =  0,  Sk+i  =  ±  [(I  —  B)  +  S2k] ,  k  =  0, 1, . . .  (9) 

with  (7  =  2—  Soo.  It  is  valid  for  x  such  that 
||  I— B(x)\\  <  1  thus  providing  an  estimate  of  the  neigh¬ 
borhood  U. 

The  relative  input-output  influence  of  each  state  in 
the  balanced  realization  is  measured  by  the  singular 


value  functions  (SVF)  (see  [2]),  which  are  the  point- 
dependent  analogs  of  the  Hankel  singular  values  for 
a  linear  system.  Figure  4  shows  the  resulting  SVFs 
for  the  pendulum  system  (3-4).  We  see  that  in  the 
balanced  realization,  one  state  has  roughly  twice  the 
input-output  influence  as  the  other.  The  values  of  the 
singular  value  functions  at  the  origin  are  0.0376  and 
0.0193,  respectively.  As  expected,  these  are  in  close 
agreement  with  the  Hankel  singular  values  of  the  lin¬ 
earized  system,  which  are  0.0384  and  0.0217,  respec¬ 
tively. 


Figure  4:  Singular  value  functions  for  pendulum  system 
in  a  small  neighborhood  of  the  origin  (denoted 
SVF  1  and  SVF  2). 


3  Conclusion 

We  have  illustrated  a  computational  approach  for  prac¬ 
tical  implementation  of  Scherpen’s  nonlinear  balancing 
procedure  via  application  to  the  forced  damped  pen¬ 
dulum  system.  A  Monte-Carlo  approach  provides  an 
estimate  of  the  controllability  energy  function.  The  ob¬ 
servability  energy  function  is  computed  via  numerical 
integration  and  is  also  suitable  as  an  empirical  method. 
Numerical  implementation  of  the  Morse-Palais  lemma 
is  achieved  using  a  successive  approximation  algorithm 
for  an  operator  square  root  which  appears  in  Palais’ 
proof.  The  corresponding  convergence  criterion  gives 
an  estimate  of  the  neighborhood  in  which  the  given 
function  is  locally  quadratic.  The  balanced  realization 
for  the  pendulum  system  results  in  one  state  having 
roughly  double  the  input-output  influence  as  the  other 
state. 
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